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ABSTRACT 

Intervening H I 21-cm absorption systems at 2; > 1.0 are very rare and only 4 
confirmed detections have been reported in the hterature. Despite their scarcity, 
they provide interesting and unique insights into the physical conditions in the 
interstellar medium of high- 2; galaxies. Moreover, they can provide independent 
constraints on the variation of fundamental constants. We report 3 new detec- 
tions based on our ongoing Giant Metrewave Radio Telescope (GMRT) survey 
for 21-cm absorbers at 1.10 < 2;abs < 1-45 from candidate damped Lyman-a 
systems. The 21-cm lines are narrow for the Zabs= 1.3710 system towards SDSS 
J0108-0037 and z^^,= 1.1726 system toward SDSS J2358-1020. Based on line 
full- width at half maximum, the kinetic temperatures are < 5200 K and < 800 K, 
respectively. The 21-cm absorption profile of the third system, Zabs=l-1908 sys- 
tem towards SDSS J0804+3012, is shallow, broad and complex, extending up to 
100 km s~^. The centroids of the 21-cm lines are found to be shifted with respect 
to the corresponding centroids of the metal lines derived from SDSS spectra. 
This may mean that the 21-cm absorption is not associated with the strongest 
metal line component. 

Subject headings: cosmology: quasars: absorption lines-radio lines:galaxies 

1. Introduction 

The diffuse interstellar medium exhibits a wide range of physical conditions such as 
temperature, density and radiation field that are infiuenced by in-situ star-formation, cosmic 
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ray energy density, photoelectric heating by dust as well as mechanical energy input from 
both impulsive disturbances such as supernova explosions and steady injection of energy in 
the form of stellar winds. Therefore, understanding the physical conditions in the gas and the 
processes that maintain these is important for understanding galaxies and their evolution. 
The damped Lyman-a systems (DLAs), with log A^(H l)>20.3, are a major reservoir of H I 
at high z and possibly the progenitors of present-day galaxies (see Wolfe et al. 2005). At 
high 2;, despite many attempts, only a handful of DLA galaxies have been detected based on 
line and/or continuum emission (see M0ller et al. 2004). 

Our understanding of physical conditions in DLAs at z > 1.8 is based primarily on the 
absorption lines of H2 molecules and atomic fine-structure lines. A systematic search for H2 
in DLAs at {z^^a^ > 1.8) has resulted in a detection in ~15% of the cases (see Ledoux et al. 
2003). Usually, the DLAs with H2 absorption also show absorption lines of C I, C I*, C I** 
and C II*. Detailed investigations show that H2 components have properties similar to that 
of the cold neutral medium (CNM) in a radiation field of moderate intensity originating from 
local star formation activity (Srianand et al. 2005). Using C ll* absorption lines Wolfe et al. 
(2004) have concluded that most of the DLAs at high z consist of CNM gas. However, above 
techniques cannot be used to probe nature of the gas at low and intermediate redshifts. In 
this case detecting 21-cm absorption line is one complementary way to probe the nature of 
absorbing gas. There have been several systematic searches for 21-cm absorption in DLAs 
undertaken by various groups (Briggs & Wolfe 1983; Lane 2000; Kanekar & Chengalur 2003, 
Curran et al. 2005) with limited success. To date, in the literature 38 DLAs have been 
searched for 21-cm absorption, resulting in 17 detections most of which occur z < 1 (see 
Fig. [1]). The column density of H I for an optically thin cloud that covers a fraction fc of 
the background source is related to the optical depth t{v) in a velocity interval v and v+dv 
and to the spin temperature (Tg) by 



The low detection rate at high z can therefore be attributed either to the gas being warm 
(high Ts) or to a low value of the covering factor (fc) through high z geometry effects (see 
Kanekar & Chengalur, 2003 and Curran & Webb, 2006, respectively). However, the redshift 
coverage is sparse and measurements are available only for a few systems at 1 < z < 2 
(Fig. [1]). To improve the statistics we have started a systematic survey of 21-cm absorption 
at 1.10< z <1.45 using the 610-MHz receiver at GMRT. Here, we report the results from 
the first phase of this survey that has resulted in 3 new detections. 

Obtaining a new DLA sample at z < 1.6 is virtually impossible in the absence of 
UV-spectrographs in space as the redshifted wavelength of Lyman-a line falls below the at- 
mospheric cutoff. However, Rao & Turnshek (2000) have shown that the DLAs can be prese- 
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lected on the basis of equivalent widths of Mg ll, Fe ll and Mg I absorption hnes. Specifically, 
they found that 50% of the absorbers with rest frame W(Fe llA2600) and W(Mg llA2796) 
greater than 0.5 A were confirmed DLAs. The detection rate becomes 100% when W(Mg l)> 
0.5 A. It is also clear from Fig. 3.4 of Lane (2000) that such a selection will also ensure that 
50% of these candidates are detected in 21-cm absorption. Motivated by this, we have be- 
gun a GMRT 21-cm survey of DLA candidates selected on the basis of W(Mg ll) in the 
redshift range 1.10< z < 1.45. Our complete sample is drawn from the catalog of 7421 
Mg II systems with W(Mg ll)> lA and 0.3< z <1.9, detected along the line of sight toward 
45,023 QSOs in SDSS DR3 (Prochter, Prochaska & Buries, 2006). There are 2857 systems 
at 1.10< z < 1.45. Out of these, we have selected all the 26 systems that have an estimated 
610 MHz flux density in excess of 100 mJy from NVSS and First catalog. These form our 
main sample. In addition to this we have also observed three other sources [J0108— 0037 
(York et al. 2006), J0240-2309 and J1604-0019 (Lanzetta, Wolfe & Turnshek, 1987)] that 
have total flux density at 610 MHz in excess of 1 Jy and have strong Mg ll systems. 



2. Observations and Results 



In the first phase of our GMRT survey we have observed 10 systems (Tabled]) from our 
sample. Usually, a 1 MHz bandwidth split into 128 frequency channels was used to acquire 
data in the two circular polarization channels RR and LL. Broader (2 MHz for J0804-I-3012) 
or narrower (0.5 MHz for J2358— 1020) bandwidths were used in the subsequent observing 
runs that were carried out for confirming the detected absorption. The data were reduced 
using the Astronomical Image Processing System (NRAO AIPS package) following standard 
procedures as described in Gupta et al. (2006). Final spectra for the three new detections 
and six non-detections are presented in Figs. [2] and [3] respectively. Analysis of Zi,hs= 1.3647 
towards J0240— 2309 is presented in Srianand, Gupta & Petitjean (2006). In most of the 
cases the background quasar is unresolved for the typical 5"x5" synthesized beam achieved 
in our observations. However, for J0804+3012, J1411-0300 and J1604-0019 the radio 
sources are extended (see Fig. H]). For J0804-(-3012, 21-cm absorption is detected in the 
spectrum extracted towards the radio peak (Fig. [2]). Higher spatial resolution as well as S/N 
observations will be required to investigate its variation across the source. In the case of 
J1411— 0300, the radio peak is not consistent with the location of the optical source. 21-cm 
absorption is not detected towards either the strongest radio peak (PI) or another peak (P2) 
northeast to PI (Figs. [3] and H]). No 21-cm absorption is detected towards any of the three 
prominent radio peaks seen in J1604— 0019. Fig. [3] shows the spectrum towards the central 
component. A narrow absorption is present with Tpeak ~0.04±0.01 and FWHM=6.1±2.3 
km s~^ near the expected frequency towards J0845+4257 (see Fig. [3]). The feature is present 
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in both the polarizations and in the spectra extracted using different basehnes and time 
ranges. At this stage we consider this as a tentative detection. 

2.1. Zabs=1.3710 system towards J0108-0037 

The 21-cm absorption profile is well fitted with a single Gaussian having Xpeak = 0.068±0.002 
and FWHM = 15.4±0.6 km s~^. From the observed optical depth we get the integrated col- 
umn density, A^(H l) = 2. Ix 10^^ {Tg / fc) cm~^. High frequency VLBA observations suggest 
structures extending over 60 mas (~ 50 pc; Ho=71 km s~^ Mpc~\ fim=0.27, Q\=0.73) with 
significant flux density in extended emission (Beasley et al. 2002). Thus fc is uncertain. As- 
suming thermal broadening, gives a gas kinetic temperature = 21.86 x (|^^)^ < 5200 
K. This is an upper limit as various other effects also could contribute to the line width. 
Now we can assume Tg = Tk as is the case for thermodynamic equilibrium. If the gas is a 
typical CNM (say = Tk = 100 K) then log [fcN(B. l)] = 20.30. Thus, A^(H l) in the 21-cm 
component is consistent with it being a DLA. The 21-cm centroid is within 4 km s^^ to that 
expected from the metal lines detected in the SDSS spectrum. The Zn ll+Cr ll blend at A 
= 2062 A, the Si llAlSOS and the Mn ll lines are clearly detected even in the low-dispersion 
SDSS spectrum. 

2.2. Zabs=l-1908 system towards J0804-K3012 

21-cm absorption is clearly detected towards J0804+3012 which is partially resolved in 
our observations (Fig. Hj). The absorption profile is consistent with a main component (with 
FWHM ~ 44 ± 19 km s^^and Tpcak = 0.006 ± 0.001) and a broad wing (extending up to 
~ 100 km s~^) towards higher frequency. Our two-epoch observations with different spectral 
resolution confirm the main component which also shows the expected Doppler shift. How- 
ever, deeper observations are required to ascertain the strength of the wing. It is interesting 
to note that the expected position of 21-cm absorption based on the SDSS spectrum coin- 
cides with the wing of the 21-cm absorption and not the stronger main component. Based 
on the profile obtained on 15—16 September 2006, we get an integrated column density of 
iV(H l)=6.5xl0^^ (Ts/fc) cm~^. Unlike the other two 21-cm absorbers, we do not detect 
Si IIA1808 (with W<0.3A) and Zn llA2026 lines (with W<0.2A) in the SDSS spectrum. 
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2.3. Zabs=1.1726 system towards J2358-1020 

Relatively weak 21-cm absorption is clearly detected and the profile can be well fit- 
ted with a single Gaussian component with rpeak= 0.035±0.004 and FWHM = 6.0±0.8 
km s-^(Fig[2]). This implies <800 K and log [/^ A^(H l)]< 20.5 in the 21-cm absorbing 
component. If the absorbing gas happens to be a CNM then log [fc A^(H l)] = 19.6. The value 
of fc could be close to unity as VLBA observations at 5 GHz by Fomalont et al. (2000) show 
the source to be < 0.5 mas, which corresponds to a projected size of < 4 pc at 2;abs=l-1726. 
Thus the inferred range in A^(H l) is consistent with the 21-cm absorbing gas being a sub- 
DLA. Like the previous system, the 21-cm absorption is shifted by 70 km s~^ with respect 
to the optical redshift based on the SDSS data. The Zn ll+Cr ll blend at A = 2062 A is 
clearly detected in the SDSS spectrum. The compactness of the background radio source 
and the simplicity of the 21-cm absorption profile make this system an ideal case for probing 
the variation of fundamental constants. 



3. Discussion 

In this letter we have reported the detection of 21-cm absorption in three systems. 
This has substantially increased the number of 21-cm absorption systems at z > 1.0, as 
prior to this only four systems were known. In two of our systems {zs,hs= 1.1908 towards 
J0804+3012 and Zabs= 1.1726 towards J2358-1020) the centroid of the strongest 21-cm 
absorption component is shifted with respect to the centroid of metal lines derived from SDSS 
spectra. This may mean that the strongest 21-cm absorption component is not associated 
with the strongest metal absorption component. Such shifts have already been observed in 
other systems (see Srianand et al. 2005; Kanekar et al. 2006; Heinmiiller et al. 2006). It is 
also known that H2 carrying components in DLAs, that are believed to trace the CNM gas, 
often show similar shifts (although less pronounced) with respect to the strong metal line 
components (Ledoux et al. 2003. H2 components often distinguish themselves with respect 
to the rest of the velocity components by the presence of excited fine-structure lines of C I 
and C II and in some cases by large depletion factors (e.g. Rodriguez et al. 2006). It would 
be very important to obtain high resolution optical spectra of our sources to see if this is 
also true in the case of the 21-cm absorption systems. 

We also report upper limits on 7 systems with one of them being a tentative detection. 
We find a detection rate of 21-cm absorption of roughly 30% in our survey so far. This is 
much higher than that achieved earlier (see for example Lane, 2000), mainly because we 
have restricted ourselves to the stronger Mg ll systems. Since the objects observed during 
the first-phase of the survey are drawn randomly from the parent sample, most probably the 
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detection rate of our whole survey will be close to what we have achieved so far. In Fig. [5] we 
compare the properties of our systems with those of DLAs in the sample of Rao & Turnshek 
(2000). The distribution of different equivalent width ratios together with the distribution 
of W(Mg ll) in our sample is consistent with what is observed in the DLA population. 

It is important to obtain independent constraints on the time variation of fundamental 
constants in order to resolve the controversy regarding the variation of a (Murphy et al. 
2003; Srianand et al. 2004). As the energy of the 21-cm transition depends on the electron- 
to-proton mass ratio (/x), the fine-structure constant [a) and the proton G-factor (Gp), high 
resolution optical spectra in conjunction with high resolution 21-cm spectra can be used to 
probe the variation of these constants (Tzanavaris et al. 2005). However, it is important 
first to understand the origin of the relative shifts that we observe between the redshifts of 
the atomic heavy element lines and the 21-cm absorption. This needs detailed modeling of 
the absorption systems taking into account all transitions simultaneously. This is what we 
plan to do in the near future. 

RS and PPJ gratefully acknowledge support from the Indo-French Centre for the Pro- 
motion of Advanced Research (IFCPAR) under contract No. 3004-3. We thank the referee 
for useful comments, Rajaram Nityananda for encouragement and GMRT staff for their 
co-operation during observations. The GMRT facility is run by NCRA-TIFR. 
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Table 1: Observational log and results. 



Name 






Time 
(hr) 


Peak 


rms** 




t" 


J0108-0037 


L 


.3710 


4.4 


1276 


2.9 


0. 


,070 


J0214+1405 


1, 


.4463 


4.4 


220 


2.7 


<0. 


,012 


J0240-2309 


1 


.3647 


8.0 


5100 


5.2 


<0. 


,001 


J0748+3006 


1, 


.4470 


5.5 


347 


3.5 


<0. 


,010 


J0804+3012 


1, 


.1908 


6.0 


2050 


1.9 


0, 


,006 








6.2^^ 


2087 


1.6 


0. 


,006 


J0845+4257 


1, 


.1147 


5.0 


224 


3.0 


<0. 


,013 


J1017+5356 


1, 


.3055 


4.4 


127 


2.3 


<0. 


,018 


J1411-0300^ 


1, 


.4160 


3.9 


244 


3.0 


<0. 


,012 










161 


3.0 


<0. 


,019 


J1604-0019 


1, 


.3245 


3.3<^-^ 


375 


3.9 


<0. 


,010 


J2358-1020 


1, 


.1726 


3.9'^-'^ 


443 


1.4 


0. 


,033 








5.5^ 


420 


2.5 


0. 


,035 



" peak flux (mJy/beam); spectral rms (mJy/beam/channel) 
peak optical depth or 1-cr limit per channel; 
2 MHz band width. 

The two spectra correspond to peaks Pi and P2. 
^ 256 channels used; ^ 0.5 MHz bandwidth 
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Fig. 1. — Line-histogram is the redshift distribution of the DLAs for which redshifted 21- 
cm observations are reported (Table. 1 of Curran et al. 2005 and z — 2.347 towards PKS 
0438—436 from Kanekar et al. 2006). Shaded histogram is for the confirmed 21-cm absorp- 
tion systems. The redshift range covered in our GMRT survey is also marked. 
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Fig. 2. — GMRT H I spectra of the sources with 21-cm absorption. Single Gaussian fits are 
overplotted in the case of J0108— 0037 and J2358— 1020. Arrows mark the expected positions 
of 21-cm absorption based on the metal absorption lines. 
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Fig. 3. — Spectra for the sources with non-detection in 21-cm absorption. Arrows mark the 
expected position of 21-cm absorption based on the metal absorption hnes. 
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Fig. 4. — GMRT maps of sources resolved in our observations with contour levels as 
nx (-1,1,2,4,...) mJy/beam. (a) J0804+3012 with an rms of 0.84 mJy/beam and beam 
of 6.36"x4.09" with position angle (PA)= 67° and n=4. (b) J1411-0300 with an rms of 
1.25 mJy/beam and beam of 5.67"x4.32" with PA=28° and n=5. (c) SDSS J1604-0019 
with an rms of 1.50 mJy/beam and beam of 7.55" x4. 34" with PA=32° and n=8. The 
position of the optical source in each image is marked as a cross. 
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Fig. 5. — Distributions of W(Mg ll) and different line ratios. Solid and dotted histograms 
are respectively for our GMRT sample and the sample of DLAs (downscaled by factor 2) 
studied by Rao & Turnshek (2000). Light and dark histograms are for the systems observed 
for 21-cm absorption so far and the systems with detections respectively. 



